Search for the decay K\ 37 
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We performed a search for the decay K\ — > 3j with the E391a detector at KEK. In the data 
accumulated in 2005, no event was observed in the signal region. Based on the assumption of K^ — > 
37 proceeding via parity- violation, we obtained the single event sensitivity to be (3.23±0.14) x 10 -8 , 
and set an upper limit on the branching ratio to be 7.4 x 10 -8 at the 90% confidence level. This is 
a factor of 3.2 improvement compared to the previous results. The results of K\ — >• 37 proceeding 
via parity-conservation were also presented in this paper. 



PACS numbers: 13.25.Es,11.30.Er,12.15.-y 

We report the first results of a search for the de- 
cay K\ — >> 37 since the last experimental update in 
1995 [1]. Although the decay is forbidden by charge- 
conjugation invariance, it can proceed via weak parity- 
violating interactions without violating CP. But due to 
further suppressions by the gauge invariance and Bose 
statistics [2|, the branching ratio (BR) of if£ — » 37 is 
expected to be very small. Assuming the decay proceeds 
via K\ — )> 7r°7r°7 — >• 37 process with two 7r 0, s internally 
converting to photons, the calculated BR(if£ — )> 37) is 
3 x 10 -19 [3] . Recently, a new calculation based on parity- 
violating model showed that the BR should be in the 
range of 7 x lO" 17 < BR(if£ -> 37) < 1 x 10~ 14 Q. 

The E391a experiment [H|, Q was conducted at KEK 
using neutral kaons produced by 12 GeV protons inci- 
dent on a 0.8-cm-diameter and 6-cm-long platinum tar- 
get. The proton intensity was typically 2 x 10 12 per spill 
coming every 4 sec. The neutral beam [7], with a solid 
angle of 12.6 /istr, was defined by a series of six sets 



of collimators and a pair of sweeping magnets aligned 
at a production angle of 4 degrees. A 7-cm-thick lead 
block and a 30-cm-thick beryllium block were placed be- 
tween the first and second collimators to reduce beam 
photons and neutrons. The beam size at the entrance of 
the detector (11.8 m downstream of the target), which 
was measured with the E391a detector by reconstructing 
the Kl H> 3tt° decay, was 3.7 cm (FWHM) including the 
effects of detector resolution. The beam line was kept 
in vacuum at 1 Pa after 5 m downstream of the target 
and 1 x 10 -5 Pa inside the fiducial decay region. The 
if £ momentum measured at the entrance of the detector 
peaked around 2 GeV/c. 

Figure Q] shows a cross-sectional view of the E391a de- 
tector and defines the origin of the coordinating system. 
The detector components were cylindrically assembled 
along the beam axis. The electromagnetic calorimeter, 
labeled "Csl", measured the energy and position of the 
photons from if£ and 7r° decays. It consisted of 496 
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FIG. 1: Schematic cross-sectional view of the E391a detector. 
"Om" in the scale corresponds to the entrance of the front 
barrel (FB) detector. K^s entered from the left side. 



blocks of 7 x 7 x 30 cm 3 undoped Csl crystal and 80 
specially shaped Csl blocks used in the peripheral re- 
gion, covering a circular area with a 95 cm radius. In 
order to allow beam particles to pass through, there was 
a 12 x 12 cm 2 beam hole located at the center of the 
calorimeter. The main barrel (MB) [8| and front bar- 
rel (FB) counters consisted of alternating layers of lead 
and scintillator sheets with total thicknesses of 13.5 Xq 
and 17.5 Xo, respectively. To identify charged particles 
entering the calorimeter, an array of plastic scintillation 
counters (CV) with a 12 x 12 cm 2 beam hole at the cen- 
ter was placed 50 cm upstream of the calorimeter. Multi- 
ple collar-shaped photon counters (CC00, CC02-07) were 
placed along the beam axis to detect particles escaping 
in the beam direction. The CC02 was located at the up- 
stream end of the K\ decay region. The CC03 filled the 
volume between the beam hole and the innermost lay- 
ers of the Csl blocks in the calorimeter. The vacuum 
region was separated by a thin multi-layer film ("mem- 
brane" ) between the beam and detector regions. Detailed 
descriptions of the E391a detector are given in 0, 

In this analysis, we used the data taken in the periods 
Run-II (Feb.-Apr. 2005) and Run-III (Oct.-Dec. 2005) 
of E391a. A hardware-based trigger system was used for 
data-taking, which required two or more shower clusters 
in the Csl calorimeter with a cluster energy larger than 60 
MeV. We imposed online cuts on the CV and some other 
photon counters. The K\ decays were simulated using 
the GEANT3 Monte Carlo (MC) generator [10] and were 
overlaid with accidental hits taken with a target-monitor 
trigger. Since the decay of K\ — > 37 is via weak in- 
teractions, parity conservation is not guaranteed. There 
were three different models considered for the simula- 



tions: the phase space, the parity- violating 11], and the 
parity-conserving [12| interactions. 

Candidates of K\ — » 37 were selected by requiring 
three photon-like clusters in the Csl calorimeter with- 
out any in-time hits in the other detectors. All clusters 
were required to be between 25 cm to 88 cm from the 
center of the beamline. An additional selection crite- 
rion on the transverse momentum of K\ — >• 37 candi- 
dates, Pt < 0.05 GeV/c, was required to suppress the 
K\ — )• 7r°7r° and K\ — >> 3tt° background events with un- 



detected photons. The decay vertex of K\ candidates 
was calculated by requiring three photons to form the 
K\ mass and by constraining the vertex to lie along 
the beam axis. The MC showed that 15% of the well- 
reconstructed K\ — )• 37 events decayed before the exit of 
CC02 (z = 275 cm). To preserve acceptance, the fiducial 
decay Z- vertex (Z K o ) region was defined to be between 
200 and 550 cm. 

Most of the backgrounds to the decay K\ — ^ 37 were 
related to the detection inefficiency of photon counters 
or fusion clusters in the Csl calorimeter. A fusion clus- 
ter is defined by two or more photons which are recon- 
structed together as a single cluster. In previous E391a 
analyses, a tight energy threshold was applied to the 
MB detector to reduce the detection inefficiency. This 
caused a major signal acceptance loss due to splash- 
back and electromagnetic shower leakage from the Csl 
calorimeter to the MB. According to the MC simula- 
tions, the undetected photons in the MB mostly entered 
the upstream region, while splash-back and electromag- 
netic shower leakage tended to enter the downstream re- 
gion. Thus, a tighter energy threshold was applied to 
the upstream region of the MB to improve the detec- 
tion efficiency, and a looser threshold was applied to the 
downstream region to keep the signal acceptance. The 
particle-hit position on the MB was reconstructed using 
the TDC information measured from both ends of the 
MB counter. The calibration and the simulation of TDC 
timing were carefully treated counter by counter. The 




FIG. 2: Distributions of fs and f% (sub-figure) in the region 
of x' > 0. The values of Je and f% were normalized by the 
total sum of /e and f% in all regions, respectively. The red 
histogram shows the distribution of single-photon clusters and 
the black histogram shows the distribution of fusion clusters. 
A significant difference between the distributions in the region 
of y' > was also observed. 

fusion clusters were mainly suppressed using the neural 
network (NN) multivariate method trained by the single- 
photon and the fusion clusters selected from if£ — >> 3tt° 
MC samples. For further suppression, the cluster was 
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divided into four regions by a radial (x 1 ) and a trans- 
verse (y f ) line crossed at the center of energy of the 
cluster. The MC study showed that the energy frac- 
tions, Je = ^%Ei (i = crystals in the cluster), and the 
weighted energy fraction, fg = E^E^ x q' 2 (q' = x' or y'\ 
in these defined regions were significantly different be- 
tween the single-photon and the fusion clusters as shown 
in Fig. El By applying cuts on and NN fusion to- 

gether, the dominant background source, fusion-related 
K\ — >> 37T° events, were completely rejected in the MC. 

Although the reconstructed Z^o of the K\ — >> 7r°7r° 
events with undetected photons was not precisely mea- 
surable, MC studies on all types of K\ — >• 7r°7r° back- 
ground events showed that the difference of the measured 
Z K o and the true Z K o had a mean of only 20 cm and 
a deviation of 10 cm. Since this difference is small, the 
measured Z x o was used to reconstruct the invariant mass 
of the i th and j th photons, m^-, where the photons were 
sorted by carried energies. Events were then rejected if 
the reconstructed rriij matched the mass of pion, m^o. 
There were three possible combinations to form the rriij , 
and the relatively significant ra^o peaks were observed in 
the m23 and mis distributions. After rejecting the events 
with the values of myi matching ra^o, the mis — ^23 dis- 
tribution of the MC is shown in Fig. [3l Events inside 
the cross region were further rejected, and then the four 
corner regions were defined as the signal region. 
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events with the NN output larger than 0.7, all the remain- 
ing 68 K% -> 3tt° MC events were rejected and 62.9% of 
37 MC events remained. The K\ — >• 77 de- 
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FIG. 3: Distribution of mi3-m23 of the MC events with three 
clusters in the Csl calorimeter. The contour shows the K\ — > 
37 parity- violating results, and the dots show the K\ — > 7v°7r° 
results. The events inside the cross region were rejected. 

In the MC study, the K\ 3tt° and K\ 37 events 
showed different behaviors in the distributions of clus- 
ter energy and position. This was due to the fact that 
the number of photons appearing in the final state of the 
two processes were different. These variables, which re- 
lied on the measurements from the Csl calorimeter, and 
Z K o were combined together by using the NN method, 
and the results are shown in Fig. 2J After requiring the 
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FIG. 4: Distribution of Cluster radial position vs. cluster 
energy of the MC events and the distributions of the NN 
output (sub-figure). The banded contour shows the if£ — » 37 
parity- violating results, and the line contour shows the K\ — »• 
3tt° results. The significant differences in the distributions 
were also observed in the other two clusters. 

cay could contribute to the background if an accidental 
cluster arrived in the calorimeter in-time with the event. 
Since the energies of accidental clusters were relatively 
small compared to decayed photons, the K\ — >> 77 was 
easily identified by requiring two higher energy clusters. 
The center of energy of the two decayed photons should 
distribute around the beam center, with combined Pt 
equal to that of their parent K^. Thus, the event was 
rejected if the center of energy of the two highest energy 
clusters was less than 4 cm from the beam center. 

With all cuts applied, 3 events of 4 x 10 9 generated 
K\ — >• 7r°7r° MC events remained in the signal region as 
shown in Fig. [3j corresponding to 0.16±0.10 events af- 
ter normalization. We generated 3.2 x 10 10 K\ — >• 3tt° 
MC events (104% of data) for Run-II and 7 x 10 10 MC 
events (322% of data) for Run-III. No events passed the 
K\ -» 37 cuts. For the K\ -» 77 MC, only one in 
4 x 10 9 generated events passed the cuts, corresponding 
to 0.03^o'o3 even ts after normalization. For the MC with 
no event passing the cuts, the number of the survived 
events was conservatively set to be 0^J. The total num- 
ber of expected background events from the three sources 
was then estimated to be 0.191q' 9 o, where the quoted er- 
ror includes statistical and systematic uncertainties. The 
total error is dominated by the K\ — >• 3tt° statistical 
uncertainty. The MC events were normalized using the 
number of data events in the region of 0.13 GeV/c 2 < 
77112 < 0.17 GeV/c 2 (normalization region). The mis 
distributions of the data and the MC results in this re- 
gion are shown in Fig. [5] The region of 0.13 GeV/c 2 < 
mis < 0.18 GeV/c 2 and mi 2 > 0.18 GeV/c 2 was defined 



TABLE I: Summary of the estimated numbers of the back- 
grounds. The quoted errors include statistical and systematic 
uncertainties. 
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FIG. 5: Distributions of mi3 of the events in the normaliza- 
tion region (0.13 GeV/c 2 < m 12 < 0.17 GeV/c 2 ). The points 
with error bars show the data, the star sign shows the only 
survived K\ — > 3tt° MC event, and the shaded histograms 
show the K\ — > 77 MC events. The hollow solid histogram 
is the sum of the K° L 3tt°, K° l -+ tv 7v° and K° L 77 MC 
results. 



as the upper sideband. The region of 0.13 GeV/c 2 < mis 
< 0.18 GeV/c 2 and m 12 < 0.12 GeV/c 2 was defined as 
the lower sideband. In the full data, we observed 164 
events with a MC prediction of 158.9 ± 8.2 events in the 
upper sideband region and 6 events with a prediction of 
7.1 ± 1.2 events in the lower sideband region. Since the 
results of the three background sources well described the 
data in both normalization region and sidebands, other 
background sources, such as neutron interactions, were 
neglected. The estimated backgrounds are summarized 
in Table H 

With all selection cuts applied to the data, no events 
survived in the signal region (Fig. [6]). The single event 
sensitivity for K\ — >• 37 was defined as 

SES(K° L 3 7 ) 



A(K° L -> 37) • N(K° L decays) ; 



where A(K^ — >• 37) is the acceptance for K\ — >• 37 and 
N(Kl decays) is the integrated K\ flux. The K\ flux 
was evaluated by the K\ — >• 7r°7r° mode and was cross- 
checked by the K° L -> 3tt° mode. The K Q L fluxes at 10 m 
from the target were determined to be (1.57±0.09) x 10 11 
for Run-II and (1.11 ± 0.07) x 10 11 for Run-III based on 
the number of decays downstream of that point. The 



quoted error in the K\ flux combines statistical and sys- 
tematic uncertainties. Systematic uncertainties from dis- 
agreements between data and the MC simulation domi- 
nated the total error. The K\ — >• 37 acceptance varied 
with decay models. Results from three models are sum- 
marized in Table [Til phase space, parity- violating and 
parity-conserving interactions. The upper limits at the 
90% confidence level were calculated based on Poisson 
statistics. The parity- violating model (CP conserved) 
was used to obtain the final result and set an upper limit 
of the BR(K£ 37) to be 7.4 x 10~ 8 at the 90% confi- 
dence level. This is a factor of 3.2 improvement over the 
previous results 
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FIG. 6: Distribution of m\z-vn<2?> of the data with all selection 
cuts imposed. No event was observed in the signal region. 
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TABLE II: Summary of the acceptances of the K\ — > 37 decay {A(K\ — > 37)), combined single event sensitivities 
(SES combined), and the upper limits (UL) at the 90% confidence level for different K° L —> 3j decay models. 



Decay Model A(K° L -> fry) SES com bined UL 

Phase space (0.99 ± 0.01) x 1(T 4 (3.75 ±0.16) x 10~* 8.62 x 10 _ * 

Parity violation (1.15 =b 0.02) x 10" 4 (3.23 ± 0.14) x 10" 8 7.42 x 10" 8 

Parity conservation (1.11 ± 0.02) x 10" 4 (3.28 ± 0.14) x 10" 8 7.54 x 10" 8 
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